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o, 

^^ , Abstract When a satellite galaxy falls into a massive dark matter halo, it suffers the dy- 

namical friction force which drag it into the halo center and finally it merger with the 
central galaxy. The time interval between entry and merger is called as the dynamical 
friction timescale (Tdf). Many studies have been dedicated to derive T^f using analyti- 
O ' cal models or N-body simulations. These studies have obtained qualitative agreements 

*— > , on how Tdf depends on the orbit parameters, and mass ratio between satellite and host 

^ ' halo. However, there are still disagreements on the accurate form of Tdf. In this paper, we 

present a semi-analytical model to predict Tdf and we focus on interpreting the discrep- 
ancies among different studies. We find that the treatment of mass loss from satellite by 
s^ , tidal stripping dominates the behavior of Tdf. We also identify other model parameters 

J^ ■ which affect the predicted Tdf. 

Key words: methods: analytical — methods: numerical — galaxies: haloes — galaxies: 
CSJ ' evolution — galaxies: interactions — cosmology: dark matter 

o 

(3 ; 1 INTRODUCTION 

In the standard cold dark matter (CDM) model, structure (dark matter halo) grows in a hierarchical man- 
ner. During the merger of two dark matter haloes, the less massive one becomes the satellitqj (or subhalo) 

k^ , of the more massive one ( host halo). The satellite will orbit in the host halo and finally merger with the 

%^ • host halo. Halo mergers play an important role in the formation and evolution of galaxies, as they can 

. . ! significantly affect the star formati on rate, colors and morphology of galaxies (e.g., lBenson et al.luCKM 

12004 ; Kang et al.. 2005; Kazantzidis et al.ll2008i) . Therefore, one inevitable question about galaxy for- 
mation and evolution in the CDM scenario is to find out how long it takes for the satellite to merge with 
the host halo. 

Dynamical friction is the primary mechanism which decreases the orbital energy and angular mo- 
mentum of satellite, and drag it to the host halo center. Description of dynamical friction was firstly 
given by IChandrasekhar I (1 1943b . who derived a formula of dynamical friction based on the idealized 
case that a rigid body moves through an infinite, homogeneous sea of field particles. For most cases, 
the satellite is moving in a finite host halo, and the dynamical friction timescale (Tdf ) of satellite is 
defined as the time interval between entry and merger with the host center. The simple application of 
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Chandrasekhar' s formula to dr i ve Td f for a rigid satellite is given by iBinney & Tremaind d 19871 here- 
after BT87) and lLacev & Cold (1 19931 hereafter LC93), and t hese formulas are widely used in the semi- 
analytical models for galaxy formation and evolution (e.g ., iKauffmann et al.l 19991 Cole et all |2000|; 
ISomerville & Primack[|l999l iNeistein & Weinmannl l201oi) . Earlv studv of iNavarro et al.l(ll995l) found 



that the LC93 formul a ca n accurately match their simulation results. However, the simulation results of 
ISpringel et al.l ( l200ll) and iKang et aU ( l2005h have indicated that the LC93 formula underestimates the 
merging timescale and overestimates the merger rate as LC93 is only valid for a rigid object, not for a 
living satellite in simulations. 

For a live satellite, one needs to take into account the effect of tidal force which leads to the mass loss 
from satellite and redistribution of mass inside the satellite. Deriving an analyt ical formula of Td f for a 
live satellite is nontrivial as one has to follow both the orbit and mass evolution. IColpi et al.l(ll999[ here- 
after C99) firstly questioned the conclusion of Navarro et al. (1995), and t hey found that tidal stripping 
can significantly inc rease Tdf. This con clusion was recently confirmed bv lBovlan-Kolchin et al.l ( |2008[ 
hereafter BK08) and lJiang et al.l (l2008l hereafter J08) using high resolution simulations. BK08 and J08 
both gave fitting formulas for Tdf, but with different dependence on orbit parameters. Their results differ 
by a factor up to 2 for eccentric orbits. Using semi-analytical model with the inclusion of tidal effect, 
iTaffoni et alj (l2003l hereafter T03) derived a fitting formula for Tdf . However, their results are not well 
tested by simulations. Moreover, the prediction of T03 is quantitatively inconsistent with the results of 
BK08 and JOS. 

In this paper, we use a semi-analytical model to study Tdf of satellite. Our main motivation is 
neither to get a consistent result with simulation or other models, nor to derive a reasonable Tdf , but to 
see how the model predictions are affected by various physical processes. This will tell us which process 
dominates the pre dicted Tdf , and how to int erpret the discrepancies am ong the previous studies. Our 
model is based on lTavlor & Babull ( 1200 ll) and lZentner & Bullockl ( l2003h , but with a few modifications. 
The paper is organized as follows. In Section|2]we review the previous results. We introduce our model 
in Section[3j and compare our model predictions with the previous work in Section]!] and we summarize 
and conclude briefly in Section]5] 

2 THE PREVIOUS RESULTS 

2.1 Set Up of Initial Conditions 

The first step of modelling the evolution of satellite is to set its initial conditions, including the orbit 
energy, angular momentum and initial position. The satellite is assumed to start its orbit at the virial 
radius, i?vir, of the host halo. It has an initial orbital energy equal to that of a circular orbit of radius 
rjRvii-, and the initial specific angular momentum of satellite is parameterized as j(0) = ejc, where jc 
is the specific angular momentum of the circular orbit mentioned above and e is the orbital circularity 
(note that < £ < 1). In the following, we use R„i to denote the initial mass ratio between the host and 
satellite halo, i.e., R^ = Mh(0)/Ms(0). 

2.2 The Previous Results 

Here we briefly review the previous studies on Tdf from analytical models or N-body simulations. Using 
the Chandrasekhar's formula, BT87 derived an expression of Tdf for satellite starting with circular orbit 
in an isothermal distributed host halo as 

-tdf,BT87 ~ -, — -T-i^mTdyn (i) 

where Tdyn is the dynamical time i?vii/K/ir, and In A is the Coulomb logarithm . 
Taking into account the dependence on the orbital circularity, LC93 obtained that 

,0.78 p 
rdf,LC93 = ^:^^^V Tdyn , (2) 
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Fig. 1 Dynamical friction timescale of sinkin g sate llite predicted by the form ulas of 
|Lace^& Cole ('19931 LC93), iTaffoni et al.l (l2003l T03), iBoylan-Kolchin et all (l2008l BK08) 
and lliang et al. (20081 J08). The six panels show the dependence on the orbital circularity e, 
and each panel shows the dependence on the initial mass ratio between the satellite and host 
halo. In the results of LC93, we adopt In A = ln(l + Rm), as used by T03, BK08 and J08. 
r] = 1.0 is used in all cases. 



Note that in the above two equations, the satellite is treated as a rigid object without mass loss. 
With help of N-body simulation, C99 derived Tdi for a live satellite as 



T, 



df,C99 



1.2e' 



0.4 ri^ 



fm In A 



V Tdyn 



(3) 



where /,„ refers to the remaining fraction of satellite mass due to tidal stripping. Note that C99 only 
considers minor mergers. It's difficult to use this formula as the Tdf depends on the presumed value for 

Jm- 

Using a semi-analy tical model, T03 derived their fitting formulas for T^f , and they were updated 
by [Monaco et al.l (l2007l) . Their model have incorporated the effect of tides, but they ignore this effect 
for the large satellite (with mass R^^^^ > 0.1). Here we omit the complex formula of T03. 

Using smoothed-particles hydrodynamical simulation with gas cooling and star formation in a cos- 
mological context, J08 fitted their results with T^i as: 



Tdf, 



0.9e" 



0.6 Rn 



J08 



0.855 lii(l + R,„) 



07Tdyn 



(4) 
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BK08 considered controlled N-body simulations for two halo mergers. They gave the fitting formula 
of T(if as 

pi. 3 

Tdf.BKos = 0.216ei-9^-— ^ii— -r^Tdyn . (5) 

ln(l + Rm) 

In Figure[T]we show the Tdf as function of satellite mass and orbital circularit}0 predicted by LC93, 
T03, JOS and BK08. For a full comparison with other results, we choose In A = ln(l + Rm) in the 
formula of LC93. It can be seen that all results show a clear trend that Tdf decreases with the increasing 
satellite mass, and increases with the orbital angular momentum and energy. However, the discrepancies 
among different studies are still remarkable. For example, the results of BK08 and JOS are longer than 
that of T03 and LC93. T03 agrees well with LC93 for large satellite (R^^ > 0.1), but disagrees for 
small satellite. The results of BK08 exhibits a steeper dependence on e than other results. 

3 MODELLING THE SINKING SATELLITE 

This section describes the dynamical evolution of satellite based on the model of lTavlor & Babull ( l200ll) : 
IZentner & B ullock (2003). In section lTTI we introduce the model for the mass distribution of dark matter 
halo. Then we describe the physical processes governing the orbital and mass evolution of satellite. 
These process can be independently implemented into the model, which allows us to investigate the 
effect of any specific process by tune its free parameter. 

3.1 Halo Properties 

The dark matter halo is a gravitational self-bound system. We express the size of halo in terms of its 
virial mass 71/vii and virial radius Rvh-, which is defined as the radius within whic h the mean mas s 
density of the halo is 200 times the critical density (pc) of the universe at 2; = (e.g.. lMo et al.l[l998l) . 
The Hubble constant is adopted to be Hq = 100/ikm s~^ Mpc^^ with h = 0.7 (BK08). The dynamical 
timescale can be described as 

where Vvir is the virial velocity of a halo. 

For simplicity, the dark matter halo is usually treated as a spherically symmetric system, and a 
simple formalism for the halo density pr ofile is the profile of singular isothermal sphere (hereafter, ISO 



simple rormalism tor the halo density pr ohle is the prohle 
profile), which can be described by (e.g.. lMo et al.l[l998l) 



and 

Mi< r)^^r. (8) 

As measured by N-body simulations, the halo density profile can be well described by the NFW 
profile (Navarro et al. 1997): 

[r/rs)[l + r/rsY 

with Ts the scale radius, and (5o the characteristic overdensity. From the definition of virial radius, we can 
find the characteristic overdensity that (5o = 200c'^/[3.g(c)], where c = i?vir/?'s is the halo concentration 
parameter, and g{x) — ln(l + x)~ x/{\ + x). For the NFW profile, the halo mass enclosed a radius r is 

M(< r) = M^J-^^ . (10) 



Thi'oughout this paper, we keep the orbital energy fixed as r; = 1.0 to reduce the free model parameters. 
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The halo co ncentration i s tigh tly correlated to its mass, and we use the median relation of c 
measured bv lNeto et an(l2007l) : 



c{M) = 4.67 



M 



10^^h~H'Ip 



-0.11 



M as 



(11) 



Note that there are still debates existing in the inner shape of the N FW profile (e.g., 
iFukushige & Makinoil200lllNavarro et allbOO^IStoeRbOOeHSpringel et al.ll2008l). Varying the shape 



of NFW profile or using other halo profiles [e.g., ISO profile; Hernquist porfile dHernquistl 1 1 990h l may 
derive a different Tdi ■ However, the simulation of BK08 indicated that using a different halo profile had 
a change in T^i of only 5% (see BK08 for more details). 

Except for Sec tion |4n where the ISO profile is adopted to compare the model predictions with the 
analytical results of LC93, we use the NFW profile in other studies of this paper. When the tidal effects 
are considered, the satellite halo has a NFW profile at the time of entering (t = 0), and this profile is 
subsequently modified due to tidal heating, as described in Section [j4l 

In our studies, we select the host halo mass as 10^^ AIq, which is the typical mass used to derive the 
Tdf (BK08, JOS, C99). We have also tested that the predicted T^i has a negligible effect on the host halo 
mass once the mass ratio Rm is fixed. 



3.2 Dynamical Friction 

The satellite will sink into the halo center by the dynamical friction force which is caused by the gravi- 
tational interaction between the satellite and the background 'field' particl es that make up the hos t halo 
(for a complete description, see BT87). This effect was first discussed by IChandrasekhar I (1 1943b . and 
the force generated by the filed particles is known as the Chandrasekhar dynamical friction. By assum- 
ing that the field particles follow a locally Maxwellian velocity distribution, BT87 gave the formula of 
dynamical friction as 



Mf 



-AirG^M^lnApir) 



erf (X) 



2X 



-x^ 



Vorb 



(12) 



where Worb is the orbital velocity of the satellite, and X = Voi-b/[\/2a{r)] with a{r) the local, one- 
dimens ional velocity d i spe rsion of the host halo at radius r, which can be solved from the Jeans equation 
(BT87. ICole & Laceyl [l99 6). Fo r ISO profile a{r) = Kir/%/2; for NFW profile, we use the fitting 
formula of a{r) from lZentner & Bullockl(l2003h . We choose the Coulomb logarithm In A = hi(l + R„i), 
as used by T03, J08 and BK08. 

The Equation (fT2] i was derived with the idealized assumption that the velocity distribution of the 
dark matter particles is Max wellian and iso t ropic. Although there are debates on whether this as- 
sumption is reasonab le (e.g., iManrique et al.L 120031 : IWilliams et al.l 120041: ISalvador-Sole et aU l2005l 
iBellovarv et aliuOOSl) . in this paper, we follow most authors (e.g., LC93: C99: T03: Zentener & Bullock 
2003: Fellhauer & Lin 2007: BK08) to adopt the Maxwellian and isotropic veloc ity distribution . There 
are also simulations showing that this assumption is a good approximation (e .g., ICole & Lacevl Il996t 
IShethLll996HSeto & YokovamaLll998tlKang et alil20ollHavashi et alil2003h . 



3.3 Tidal Mass Stripping 

For a live satellite, the tidal force from the host halo will strip its mass. The tidal radius, rt, is the distance 
from the center of satellite to the radius where the external differential force from the host halo exceeds 
the binding force of the satellite. The tidal radius can be simply solved from the following equation 
(Ivon Hoerneiill957HKinglll962l:lTavlor & BabulLl200l : 



w^ 



GMs{< n) 

G[2Mh{<r)/r^-47Tph{r)] ' 



(13) 
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with uj the angular speed of the satellite and ph (r) the density profile o f the host halo. The mass outside 
rt becomes unbound and is stripped gradually. iTaylor & Babull (12001 1) suggested the unbound mass to 
be stripped at the rate that 

dMs Ms{> rt) 



dt T, 



(14) 



orb 



with Torb the instantaneous orbital period (i.e., Torb = 27r/cj), which is assumed as the mass stripping 
timescale. 

There are some uncertainties in the above mechanisms of mass stripping, (i) The tidal radius cannot 
be characterized by a single radius, as the zero-velocity surface (the surface defined by the tidal radius, 
see BT87) is not spherical, (ii) The perturbation of particles within the satellite may lead to the scatter 
in UJ, and the zero-velocity surface is actually a shell of 'non-zero' thickness, while this effect is ig- 
nored in Equation ( fT3] l. So the solution of Equation ( fT3] ) is only an approximation for the tidal radius, 
(iii) The stripped mass from a satellite still remain in the vicinity of the satellite, and the interaction 
between the stripped a nd unstripped mass will perturb the satellite orbits and affect the mass loss (e.g.. 



between the stripped a n( 
lFellhauer&LinLl2007l) . 



Owing to these uncert ainties, numerical sim ulat ions have debated on h ow fast the unbound mass is 
stripped from the satellite. [Zentner et al.l ( l2005l) and iDiemand et al] ( |2007|) found a stripping timescale 
3.5 and 6 times shorter than Torb, respectiv ely. It was also po i nted out that the stripping timescale i s 
dependent on the satellite internal structures dKazantzidis et al.ll2004l : iKampakoglou & Bensonll20()7l) . 
In general, the mass loss rate can be described using a free parameter a as: 

dM, M,(> rt) 

-7—= -a — 7f , (15) 

at 1 orb 

where a describes the efficiency of tidal stripping. In Section l4T2l we will show how the T^f depends on 

a. 

3.4 Tidal Heating 

During the pericentric passage of satellite orbits, the gravitational field changes rapidly, and this in - 



uurmg tne pericentric passage or satellite orbits, tne gravitational neig changes rapiaiy, ana tnis in - 
duces a gravitational shock that can add energy to the satellite (e.g.. lGnedin & OstrikeA 1997 , 19991) 
This effect is called the tidal heating. It has been found from iV-body simulations (e.g.. lHavashi et al] 



This effect is called the tida l heating. It has been found from iV-body simulations (e.g.. lHavashi et aU , 
2003 : Kravtsov et all |2004|) that tidal heating will expand the satellite and reduce its inner mass profile. 



Havashi et alj (|2003|) introduced a modified NEW profile to describe the density distribution of a tidally 
heated satellite according to 

P^"^") "" 11/*/ — ^PNFw{r) , (16) 

where 

Ig ft = -0.007 + 0.35a:„ + 0.39a;,^„ + 0.234> , (17) 

and 

Ig — = 1.02 + 1.38a;,„ + O.STa;^ . (18) 

In Equation ( fTSI l, pnfw {r)\s the original NEW density profile of the satellite at the time of entering (< ~ 
0), ft describes the reduction in the central density of the satellite, and rte is the 'effective' tidal radius 
that describes the outer cutoff imposed by the tides. In Equation ( [TtI i and ( fTSl l. Xm = ^^[M s{t) / M s{i))] 
is the logarithm of the remaining fra ction of satellit e mass , and r^ is the scale radius of the satellite with 
NEW profile at i = 0. As shown bv lHavashiet alj (|2003|) . ft and rte are well fitted by the function of 
Xm- Both ft and rte decrease with time while a satellite is losing mass. 
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3.5 Orbital Evolution 



Here we present explicitly the equations to solve the orbit [x(r, 6)] of the satellite under gravity and the 
dynamical friction. The equation of motion for the satellite is given by 



dt2 



GMhi< r) 






(19) 



with Mh{< r) the mass of the host halo inside of radius r, and Fdf the dynamical friction force given 
by Equation ( fT2b . The orbital energy and angular momentum of the satellite will decay due to the 
dynamical friction as it is always opposite to the direction of motion. We define satellite to be merged 
with host center when it loses all its angular momentum, and Tdf is the time interval between accretion 
and mergerlj (as used also by BK08). The equation of motion and Equation dTSI l are solved using the 
fifth-order Cash-Karp Runga-Kutta method, in which an adaptive step-size control is embedded. 

4 RESULTS 

4.1 Examination on a Rigid Satellite 
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Fig. 2 The dependence of Tdf on satellite mass (left panel) and orbital circularity (right panel) 
for a rigid satellite, where Tdf is normalized to its value when Rm = 20 and e = 1, respec- 
tively. The model with a = {) means that a rigid satellite is considered. The results in red solid 
lines are computed with NEW profile, while the blue dashed lines show the results with ISO 
profile. Both model predictions match well with LC93's (balck dotted). 



Eirstly we validate our model by comparing the predicted Tdt with the LC93 result for a rigid 
satellite. LC93 derived Tdf using Equation ( fT2b and ISO profile for the host halo. In our model, we 
simply set a = to 'close' the tidal stripping and tidal heating effect, and we model the host halo with 
both NEW profile and ISO profile. 

In Eigure |2] we show the Tdf as a function of R^^^^ and e for a rigid satellite, with Tdf normalized 
to its value when R,n = 20 and e = 1, respectively. As indicated, our results in NEW (red solid) and 
ISO (blue dashed) model both have the same dependences as predicted by LC93. On the other hand, the 



^ Some (e.g.. lKravtsov et an . l2004HZentner et alil2005h define satelHte to be merged with the host halo when its distance to 
the host center is less than a fiducial radius. We find that different definitions have no significant effects. 
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Fig. 3 Comparison of Tdf between our model (a = 0) and LC93 (dotted) for a rigid satellite. 
The Tdf with NFW profile (red solid) and ISO profile (blue dashed) both agree well with 
LC93's prediction, and they agree for all orbital circularity, although only e with 0.1, 0.5, 1.0 
are given here. 



amplitudes of Tdf from the models also agree well with the results of LC93, which is demonstrated in 
Figure|3] The difference resulted by varying halo profile are small and negilible, which is also concluded 
by BK08. 



4.2 Dependence on Tidal Stripping Efficiency a 

In this section we study the effects of tidal stripping efficiency (a). In Figure|4]we show the predicted 
Tdf with different values of a. A larger value of a corresponds to a stronger tidal field or a rapid mass 
loss from the satellite. The results show a remarkable trend that the Tdf is increased when the tidal field 
becomes stronger The reason can be seen from Figure |5] which shows the evolution of satellite mass 
and specific angular momentum with dependence on a. The initial conditions are set as Rm = 10, 
£ = 0.5, and rj — 1.0. The left panel shows that a stronger tidal field will induce more mass loss from 
the satellite, and this effect is more distinct at the beginning. As seen from Equation (fT2] i. the amplitude 
of dynamical friction has a strong dependence on the mass of satellite (Tdf (x Af^). So a stronger tidal 
stripping will lead to a slower decay of satellite angular momentum and result in a longer dynamical 
friction timescale, as shown in the right panel of Figure |5] 

As shown in Figure [T] the predicted Tdf from the previous results disagree with each other quan- 
titatively. We believe that the main discrepancy is resulted in the treatment of tidal stripping, and we 
discuss it in more details in the following. 

- T03 ignored the tidal effects for massive satellite (with mass R^^^ > 0.1), and so their Tdf are 
consistent with LC93's. But T03 predicted a longer Tdf for low-mass satellite which suffers from 
tidal stripping. 

- The Tdf inferred by JOS and BK08 are longer than that of T03. This is because T03 adopted a tidal 
stripping efficiency that is different from those in N-body simulations. T03 al so used Equa tion ( fTSl ) 
to describe the mass loss, but with a = 1.0 which is too low. As shown bv IZentner et al.l ( l2005b . 
a higher value that a = 3.5 is required to better fit the satellite mass function from simulations 
(also see lGan et al.L 120101) . A higher value of a is also favored from Figure|6]where we compare the 
evolution of satellite specific angular momentum from our model (solid lines) with the simulation 
results of BK08 (dashed lines). We find that a = 2 can better match the simulation results. Thus the 
lower value of a used by T03 explains why they obtained a lower Tdf. 
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Fig. 4 The dynamical friction timescale (T^i ) predicted by our model. The lines with increas- 
ing thickness show the effects of tidal stripping efficiency {a ~ 1,2, 3). The Tdt is a strong 
function of a. 

- The Tdt of JOS is longer than that of BK08 for eccentric orbit (i.e., low efl The simulation of 
JOS includes the process of gas cooling and st ar formation. The halo of a sate l lite is expected to 
contract in response to the cooling of gas (e.g.. lGnedin et al.ll2004l: lAbadi et aUl2010l) . During the 
pericentric passage, the satellit e with halo contraction is resistant to the s trong tidal field, and will 
survive for a longer time (e.g.. [Weinberg et aU l2008l iDolag et al.l 120091) . Instead BK08 performs 
a higher resolution simulation, in which the satellite can avoid the artificial mass loss due to the 
numerical effects. So the satellite will deposit more mass in the eccentric orbit and suffer stronger 
dynamical friction. 



4.3 Dependence on Orbital Circularity e 

The previous results showed similar dependence of T^t on the initial satellite mass, but very different 
dependences on the orbital circularity [Equations (l2]i-(|5]l]. For example, BK08 found an exponential 
dependence of Tdf on the the orbital circularity, while others found a power-law dependence. Here we 
investigate this problem using our model with a = 1. We compute the Tdf as a function of e for a minor 
merger (Rj^/ = 0.05) and a major merger (R,^^ = 0.3), as shown in Figure]?] We find the dependence 
for the minor merger can be fitted to a power law, Tdf oc £*''*, as predicted by C99 (long-dashed). For 
the major merger, the dependence is close to the result of BK08, who found the exponential law that 
Tdf (X exp(1.9e). It is not a surprise as C99 only considers minor mergers while BK08 has more samples 



^ The results of BK08 and JOS also differ for small satellite with large e, of which the Tdf , however, are extrapolated by their 
formulas and exceed the Hubble time. 
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Fig. 5 The evolution of satellite mass and specific angular momentum (both are normalized 
to their initial value) as function of tidal stripping efficiency a. The initial conditions are that 
Rni = 10, e = 0.5, ?7 = 1.0. Strong tidal effects reduce the amplitude of dynamical friction 
and decelerate the loss of angular momentum. 
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Fig. 6 Comparison of the evolution of satellite specific angular momentum between our 
model and BK08, with three cases of initial masses and orbits as indicated. The solid lines 
with increasing thickness are the model results with a = 1,2,3, while the dashed line is the 
result from Figure 1 of BK08. The tidal stripping efficiency in the simulation of BK08 should 
be stronger than that in a model with a = 1. 



for the major mergers. Thus we argue that the dependence on orbital circularity is mainly determined 
by the distribution of mass ratio between the satellite and host halo. 



5 CONCLUSION AND DISCUSSION 

hi this paper, we study the dynamical friction timescale (Tdf ) of a sinking satellite into a host halo. 
Previous results using analytical models or simulations generally agree that the T^{ is correlated with 
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Fig. 7 The dependence of Tdf on orbital circularity for a live satellite (a = 1). For a minor 
merger (triangle), with R~^ = 0.05, the result shows a power-law dependence, which is 
similar to that of C99 (long-dashed), while for a major merger (square), with R~^ = 0.3, it 
indicates an exponential dependence, which is close to that of BK08 (short-dashed). 



the mass, orbital circularity and energy of the satellite, but disagree on the amplitude of Tdf and the 
dependence of Tdf on orbital circularity. It was unclear what contributes to these discrepancies among 
different studies. 

Aiming at inte r pretin g the se different dependences, w e use a semi-analytical model similar to that 
of iTaylor & Babull J2OOII) and IZentner & Bullockl (l2003h to derive the Tdf . Our model considers the 
main physical processes governing the evolution of satellite: dynamical friction, tidal stripping, tidal 
heating and merger. All these process are independently described by free parameters, and it allows us 
to investigate the dependence of Tdf on any process. 

Firstly, we apply our model to a rigid satellite by 'turning off the tidal stripping and tidal heating 
(i.e., a = 0). The model predictions agree well with the LC93's result on the amplitude of Tdf and its 
dependences on satellite mass and orbital circularity. Then we study the dependence of Tdf on the tidal 
stripping efficiency. We find that the Tdf depends strongly on a, with the trend that the Tdf increases 
with increasing a. A higher a leads to rapid loss of mass from satellite, than decreases the dynamical 
friction force. Thus this results in a slower decay of angular momentum and a longer Tdf . We believe 
that the main reason for the diversity of previous result is the treatment of tidal stripping. 

We also study the dependence of Tdf on the orbital circularity (e). We find that for low mass- 
ratio mergers (AIs/Mh < 0.1), Tdf is a power law of orbital circularity. While for massive mergers 
{Ms/Mfi > 0.1), the dependence of Tdf on orbital circularity is expoential. Thus we argue that the 
dependence on e obtained by different studies is determined by their samples, in which the mass ratio 
between satellite and host halo is crucial. 

In this paper, we do not model the effects of baryon, as it is difficult to include the physical processes 
governing galaxy formation, and it is still not clear how dark matter halo will respond to the baryon at 
the host halo center 

The major effect of baryon is to modify the density profile of dark matter halo. There are still debates 
about how the baryon will change the central concentration of halo. Some found that c entral density 
increa ses (e.g.. iBIumenthal et al.lll986t iGnedin et al.l |2004|) . but some disagreed with it. iGnedin et al.l 
(I2OO4I) found that the halo will become more concentrated as baryons condense in the radiative cooling. 
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and the contraction of halo is dependent on the amount of baryon. While lAbadi et al.l (1201 Ol) found that 
the response of halo contraction depends not only on h ow much baryon m ass has been deposited by 
the halo, but also on the mode of its deposition (also see lTissera et al.Ll2010l) . They sho wed that strong 
feedback by supernovae can significantly decrease the central density of halo (also see IPedrosa et all 
I2OO9I: iGovernato et al.l 120101) . The variation of T^f is about 20% when Csat/chost changes between 1 
and 2 (T03; BK08). 

There are also some studies showing that the dark matter haloes have constant density cores 



(e.g. jGilmore et al.Ll2007llde Blok et al.ll2008llKuzio de Narav et allbOOgllGebhardt & Thomasll2009l 



Hernandez & L ed I2010I). which can sig nficantly suppress the effect of dynamical friction (e.g., 
Sanchez-Salce do et al.l l2006l llnouel 1200 9'). Howerer, t he typical size of the constant density core in 
the dark matter halo is usually less than 1 kpc (e.g., Ide Blok et al.l l2008l) . The effect of the con- 
stant density core may be re markable for the evolution of globular clusters in a dwarf galaxy (e.g.. 



may be rer 
"ail 120061) . 



ISanchez-Salcedo et al.l 120061) . but not for the evolution of satellite halo in a Milky- Way sized halo. 
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